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Quantum-corrected equations of motion generically contain higher time deriva- 
tives, computed here in the setting of canonically quantized systems. The main ex- 
ample in which detailed derivations are presented is a general anharmonic oscillator, 
but conclusions can be drawn also for systems in quantum gravity and cosmology. 



Quantum effects, generically, are non-local in time, captured in effective actions or equa- 
tions by correction terms containing time derivatives of higher than second order. The path- 
integral formulation provides an intuitive explanation: entire paths connecting two points, 
CN 1 not just a local neighborhood on a single trajectory, determine observable properties. In 

gravitational theories, one expects quantum corrections with higher time derivatives for 
an independent reason: The theory is interacting and therefore should receive non-trivial 
{Sj ! quantum corrections. The only generally covariant extension of the Einstein-Hilbert action 

by functionals of the metric is by higher powers and contractions of space-time curvature 
tensors, most of which introduce time derivatives of higher than second order in equations 
of motion. 

Many properties of classical and quantum gravity are best described and analyzed in a 
canonical formulation, especially when gauge issues play a role. In canonical quantizations, 
however, it is not all too clear why and how corrections with higher time derivatives 
should result. Equations of motion with higher time derivatives imply additional degrees of 
freedom because more initial data must be provided compared to usual second-order ones. 
In a perturbative setting, the solution space does not increase in size because the surplus 
solutions are not analytic in the perturbation parameter that multiplies higher-derivative 
terms, and therefore must be discarded for self-consistency [TJ. However, terms that contain 
higher time derivatives do modify the classical solutions of second-order equations and 
therefore the new degrees of freedom they come along with play an indirect role. In 
canonical quantizations, however, one replaces the classical phase space by a set of basic 
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operators of the same number, without an obvious place for new quantum degrees of 
freedom. A clear identification of such variables together with a systematic procedure 
of deriving quantum corrections in which they appear is of general interest, not just in 
quantum-gravity research where higher-curvature corrections could be computed by such 
means. 

Effective methods for canonical quantum systems do exist, using a systematic analysis 
of quantum back-reaction of fluctuations and higher moments of a state on the evolution 
of expectation values j2j. Applying this scheme to anharmonic oscillators with potential 
V(q) = \muj 2 q 2 + U(q), it has been shown that results equivalent to those of path-integral 
based low-energy effective actions [5] are obtained. So far, these calculations, in [3] as well 
as [2], have been restricted to the first order in a semiclassical expansion by h and second 
order in an adiabatic expansion, analogous to a derivative expansion. To these orders, in 

r eS \q(t)} = / dt - m + — ^ q 2 

[ K h J \2 \ 32m 2 ^ (1 + U"{q)/mu 2 f /2 ) 

-_^V-^)- T (i + ^J J (i) 

one can see corrections by an effective quantum potential as well as a correction to the 
particle mass, but no higher time derivatives result. 

Nevertheless, the scheme provides a natural candidate for quantum degrees of freedom 
analogous to new degrees of freedom in corrections with higher-time derivatives [I] : fluctu- 
ations and higher moments of a state. As in perturbative higher-derivative theories, these 
degrees of freedom play an indirect role when an adiabatic expansion is used, because 
their equations of motion can be solved and solutions can be inserted into equations for 
expectation values to determine quantum corrections. In this article, we push the required 
expansions to higher orders to compute several new correction terms for the same systems, 
general anharmonic oscillators, and confirm that higher time derivatives appear. These 
results are collected in Sec. E] and put together in Sec. HI after our review of canonical 
effective techniques in Sec. |5J 

Looking at the details of our analysis, we will also be able to draw several general 
conclusions about properties of effective canonical dynamics. These statements, together 
with a general discussion of the relevance of our findings for (loop) quantum cosmology, 
can be found in the concluding section [5j 

2 Anharmonic oscillators 

The classical Hamiltonian for an anharmonic oscillator is given by: 

H(q, p) = t^-P 2 + \mu?q 2 + U (q) (2) 

with the particle mass m, harmonic frequency u, and an arbitrary function U(q) which, for 
many purposes such as quantum stability, is often restricted to be bounded from below. We 
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will mainly be thinking of a polynomial U (q) of higher than second order, whose total order 
is even if boundedness from below is required. (Our effective equations will be meaningful 
even when this condition is violated.) The frequency u is uniquely determined only if one 
requires U(q) to have no quadratic contribution. 

The Hamiltonian (J2]) can straightforwardly be quantized, without factor ordering am- 
biguities. Quantum states then satisfy the Schrodinger equation 

ifi^l = Hm = —p 2 \q) + -moo 2 q 2 m + U(q)\V) . (3) 
at 2m 2 

This equation takes the form of a differential equation when a representation of states, for 
instance as wave functions of q, is chosen. Effective equations, however, are independent of 
this representation choice. (There may be inequivalent representations not related unitarily, 
for instance on a non-separable Hilbert space j5]. In this case effective equations would 
depend on which representation is chosen; see the example [6] in quantum cosmology.) 
Instead of representations of wave functions, we use the general evolution equation 

djO) = ([Q,H]) 

dt ih 1 1 

for expectation values of observables O. For a Hamiltonian as given here, ()!]) applied to q 
and p gives rise to Ehrenfest's equations 

resembling the classical ones but also exhibiting quantum corrections in the force term 
— (V'(q)) compared to —V'((q)). Canonical effective equations compute the difference of 
—(V'(q)) and —V'((q)) in a systematic way. (In [7], these equations are used to prove that 
quantum mechanics has the correct classical limit for h — > 0.) 

Unless the potential is at most quadratic, Eqs. (jSJ) do not provide a closed set of equa- 
tions that could be solved for the expectation values, starting from some initial values. A 
cubic term \q 3 in the potential, for instance, gives rise to —3X(q 2 ) = — 3A((g) 2 + (Ag) 2 ). 
The second term is a quantum correction to the classical force — 3Ag 2 , but it depends on 
the position fluctuation Aq which is independent of the expectation value (q). Ehren- 
fest's equations therefore cannot be solved in this case, unless one already knows how the 
quantum state or at least its position fluctuation evolves. 

To provide a complete set of equations, effective techniques enlarge the set of Ehrenfest's 
equations by deriving differential equations for fluctuations and all moments 

G a > n :=((q-(q)Y- a (p-(P)) a ) Wcyl (6) 

where the subscript "Weyl" indicates Weyl (or totally symmetric) ordering. Since these 
variables are defined by expectation values, equations of motion for them can be derived 
from the same equation (jlj) as used earlier. Moreover, these variables provide infinitely 
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many quantum degrees of freedom independent of expectation values, just the degrees 
of freedom which, at least qualitatively, should be related to implications of higher time 
derivatives. 

Unlike expectation values, moments cannot take arbitrary values. For n — 2, they must 
satisfy the familiar uncertainty relation 



and there are analogous, but less familiar relations at higher orders. All these inequalities 
follow from the Schwarz inequality. Only if they are obeyed can the moments correspond 
to a state. If they are, the state may be pure or mixed; if a selection of pure states is 
required, additional conditions must be imposed. (See also the examples in [8].) 

Instead of calculating all the commutators required for equations of motion (j3J) of 
moments, it is usually more straightforward to take a phase-space point of view. We first 
consider the right-hand side of (jl]), which in classical mechanics would be the Poisson 
bracket of O with the Hamiltonian H. This comparison motivates the definition of a 
Poisson bracket between expectation values of arbitrary operators, 



This bracket is antisymmetric, linear and satisfies the Jacobi identity by virtue of those 
properties realized for the commutator. If we extend it to products of expectation values 
by using the Leibniz rule, we obtain a well-defined Poisson bracket on the quantum phase 
space, whose elements are states parameterized by their expectation values of q and p 
together with all moments. (The quantum phase space is symplectic. However, in effective 
equations one truncates the system at finite orders of h, which translates to finite orders of 
the moments. The Poisson tensor remains well-defined on these subspaces, but in general 
is no longer invertible. Effective equations therefore cannot be described by symplectic 
techniques.) 

Combining fl4]) and ([8]), the Schrodinger flow is described equivalently by a Hamiltonian 
flow on the quantum phase space, generated by the quantum Hamiltonian Hq({c[), (p), G a,n ) : = 
(H) (s) t (p\ ,G a ' n - The subscript indicates that the expectation value is taken in a state char- 
acterized by the values (q), (p),G a,n ; the result then defines the value of the quantum 
Hamiltonian at the quantum phase-space point corresponding to the same state. 

As a function of expectation values and moments, quantum Hamiltonians can be com- 
puted by writing (H(q,p)) = (H((q) + (q — (q)), (p) + (p — (p)))) and expanding in the 
"small" quantities q — (q) and p — (p). This formal expansion is a shortcut for a direct 
computation of the expectation value. Inserting a Taylor expansion and assuming H(q,p) 
to be Weyl ordered (in case there is any factor-ordering choice), we obtain 




(7) 



{(A), (B)}: 





1 /n\ d n H(q,p) d 
a\ \a J dp a dq n ~ a 



•a,n 



(9) 
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with the moments ([6]). (Here and from now on we identify q := (q) and p := (p) to simplify 
our notation.) 

The moments may be written in dimensionless form as G a,n = h~ n ^ 2 (moj) n ^ 2 ~ a G a ' n to 
facilitate future expansions. In terms of these dimensionless variables, we obtain, for the 
given anharmonic oscillator, the quantum Hamiltonian as: 

H Q = ^-p 2 + \muj 2 q 2 + U(q) + ^(G°< 2 + G 2 ' 2 ) + -{%/muj) n ' 2 U^ <(q)G°' n . (10) 

n=2 

The first quantum correction, depending only on fluctuations but not on expectation values, 
is a zero-point energy. The sum, on the other hand, contains products of expectation values 
and moments if there is an anharmonic potential, and therefore describes the coupling 
between quantum variables and expectation values, or quantum back-reaction. (Note that 
the moments vanish identically if n — 1. The sum therefore starts at n = 2.) 

If we look at the equations of motion for expectation values and moments generated by 
the quantum Hamiltonian [2], 

p 

Q = {q, Hq} = — 
m 



P 



oo 1 

{p, Hq} = -mu 2 q - U'(q) - V ^-(h/muj) n/2 U^ +1 \q)G ' n (11) 

^ — ' 77 .1 



n=2 



U"(a)n 

G a > n = -acoG a ~ 1 ' n + (n - a)uG a+l > n n 2 ) 

muj 

VfjaU"'{q) C a- 1 , n - lc 0,2 + h aU ""(q) C a-l,n-l C 0,3 

2(mu) 3 / 2 3!(mw) 2 

_a f VhU"'(q) 1|W+1 hU""(q) 1-raH 
2 I {mufl 2 3(mu) 2 

(not all terms are written in the last equation) we can already see that moments are 
related to higher time derivatives: Eq. (llljl can be interpreted as identifying an infinite 
linear combination of moments with the second derivative of q, in a way that also depends 
on q itself. Taking further time derivatives of the whole equation (fTTj) and inserting (TT2]) 
relates different combinations of the G a,n (no longer linear) to time derivatives of q of 
higher than second order. It is therefore clear that moments in the canonical setting play 
the role of higher time derivatives in a Lagrangian one. But so far the identification is not 
very direct, and the equations we obtain for higher time derivatives in terms of moments 
are difficult to invert. In the rest of this paper, we work out a systematic method, using 
two expansions as in [2] , to write ( fTTT) as an equation corrected by higher-derivative terms, 
eliminating the moments. 
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3 Semiclassical and Adiabatic Expansions 



Our first expansion is a semiclassical one. In a semiclassical state, the moments by def- 
inition obey the hierarchy G a,n = 0(h n ^ 2 ) so that an expansion by h to a given finite 
order makes use of only finitely many moments. Thanks to the definition of dimensionless 
variables G a,n = h~ n ^ 2 (muj) n ^ 2 ~ a G a ' n , suitable powers of % already appear as factors in 
equations of motion such as (fTTT) , and we only need to truncate the sum. 

Although the leading order of h is split off the moments when using dimensionless 
variables, each moment, as a solution of (TT2"|) . could still have higher-order corrections in 
h. For full generality, we therefore make these terms explicit by expanding by powers of 
\ffi: G a ' n = ^ e G^ n h e ^ 2 . The coefficients G® ,n are then independent of h. Some features 
of this expansion may look unfamiliar, which we explain by two comments: 

• As written explicitly, we should expect half-integer orders in % because a moment of 
order n behaves as 0{h n ^ 2 ). If all odd-order moments vanish, which is the case in a 
large subclass of semiclassical states, only integer orders appear, as naively expected. 
Such an assumption can always be made for an initial state, but non-trivial quantum 
back-reaction can easily generate non-vanishing odd-order moments. (See for instance 
the example in [9].) 

• Since H has non-trivial dimensions, the coefficients G^ n have different dimensions for 
different e. If this feature is unwanted, one can use coefficients or parameters in the 
anharmonicity potential U(q), together with m and u, to define a parameter L of 
the same dimensions as h and expand by ^h/L. The parameters of the harmonic- 
oscillator Hamiltonian have already been used to absorb the dimensions of G a,n , and 
they do not allow a combination with the dimensions of h. In the harmonic case, an 
expansion by \Jfi/ L cannot be done, and it is not necessary because the equations 
of motion for moments can then be solved exactly, showing that each moment G a ' n 
is exactly proportional to h n ^ 2 and stays so at all times. This property is no longer 
realized for anharmonic oscillators, but then there are additional parameters in the 
potential that can be used to define a suitable L. We refrain from doing so here 
because we work with a general anharmonicity U(q). Its derivatives then provide the 
correct dimensions. 

With the semiclassical expansion in H, the equations of motion for the moments partially 
decouple. At O(Ti ) and 0(h 1/2 ), we have 

G' a,n ^ya—l,n . / \ ^ya+l,n U (^JO a _ i n . . 

o = — auG + [n — a)uG G (13) 



s^a+i.n U h \q)a „ a _i n U"'(q)a ^2^0-1,71-1 



U"'(q)a f G a^ n+ i _ (a - l)(a-2) ga „3 in ,3 N . 



2(mw) 3 / 2 V U 12 
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which we could try to solve order by order. The equations for Gq ,ti couple only the n + 1 
moments at fixed order n and are linear in the moments. The equations for G"' n (and 
similarly for higher orders in h) are inhomogeneous but also contain only moments G°^ n of 
the same order n, as well as non-linear inhomogeneous terms of different orders. 

By the semiclassical expansion, the infinitely coupled original system has been reduced 
to finitely coupled subsets. In principle one could solve this system order by order, but since 
coefficients of the differential equations also depend on q for a non-trivial anharmonicity, 
they are difficult to solve explicitly. We therefore make use of a second expansion, an adi- 
abatic one, which reduces the differential equations to algebraic ones. This approximation 
will also be crucial to bring out the nature of moments as higher time derivatives. (In the 
conclusions we will comment on the nature of moments in regimes in which no adiabatic 
expansion is possible.) 

The adiabatic expansion is defined by replacing all time derivatives in equations of 
motion for moments (q and p are not assumed to change just adiabatically) by d/dt — > 
Xd/dt, expanding all coefficients Gf n = Yli^i ^e'?^ m ^> solving equations order by order 
in A, and setting A = 1 in the end. The expansion is formal because there is no guarantee 
that the series converges for A = 1. Moreover, the parameter A, unlike h in the semiclassical 
expansion, has no physical meaning. As we will see later, the procedure rather serves as a 
systematic way of organizing the appearance of derivatives of different orders. 

At zeroth order of the adiabatic approximation of G a,n = {G a,n , Hq}, we have equations 

= {G a e ^H Q } (15) 

which are algebraic rather than differential. At higher orders, 

G^ = {G^Hq] (16) 

contains time derivatives, but if we proceed order by order, we can assume that G a ^ has 
already been solved for, starting with the algebraic equation for G^q. With the time 
dependence on the left-hand side of ({TBI) known, the equation again reduces to an algebraic 
one for G a ^ +l . Proceeding order by order, the main equations to be solved are algebraic. 
(Some differential consistency conditions also arise, as we will see explicitly.) 
Combining both expansions, our moments read 

Q a,n = G a ^% e ' 2 \\ (17) 

e i 

Equations to be solved for the coefficients G^'" show both advantages noted above: They 
split into finitely coupled sets, and are algebraic. 

We now proceed to computing explicit solutions to several orders. But first we interject 
a comment on our designation of ^-orders to avoid confusion. For most of this paper, we 
will be looking at the moment equations (|12p and their solutions, and therefore speak of the 
0(fr°)-order when all fr-terms are dropped, of the 0(fr 1//2 )-order when terms linear in \fh 
are kept, and so on. These orders give us the relevant information about state properties at 
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the corresponding orders. However, if we use these moment solutions to compute correction 
terms for effective equations of expectation values, inserting the moments in (iTTj) as we 
will do in the end, the ^-orders shift because ( TTTT) contains explicit factors of H. Somewhat 
counter-intuitively, even the 0(/i°)-order in the moments will then contribute to non-trivial 
quantum corrections. This intermingling of the orders cannot be avoided because equations 
of motion for different variables — expectation values or moments of different orders n - 
have their own arrangements of ft-terms. 

3.1 Adiabatic Approximation at 0(h°) in the Moments 

At zeroth order in the adiabatic approximation, we can ignore all time dependence: {Gg'o j Hq} 
0. At zeroth order also in the \fh expansion, we have from f|T3|) 

0/~ia— l.n i / \ /nra+l.Ti U yQ)® /~ia—l t n /-i o\ 

= — auG 00 ' +(n — a)uJ& Q' ^oo U°) 



mu 



for < a < n, which gives a solution of the form 



U"(q)\ (2a 



G W- Cn ((n-a)/2)\(a/2)\ + J (19) 

with some coefficients C n if both n and a are even. Otherwise, Gq'q = 0. (For odd a, 
Eq. (1181) used with a = implies that G 'q = 0, which upon recurrence to a = 2k + 1 with 
integer k implies that Gq'q = for odd a, no matter whether n is even or odd. For n odd 
and a even in Gq,'o) Eq. (fT8|) evaluated for a = n is meaningful, with a zero value implied 
for Gq Q 1,n , only if all Gq'q with odd n and even a vanish.) 

Only the values (TT9]) with even n and a can be non-zero. Based on the zeroth-order 
equation ffT3l) alone, the C n could depend on q as well, but this possibility is ruled out by 
a consistency condition obtained at first adiabatic order [21 0] . The values of C n in general 
effective equations remain free (provided the resulting moments satisfy the uncertainty 
relation) and parameterize different choices of adiabatic states. A prominent choice is 
the anharmonic vacuum, whose values can be obtained by requiring that the harmonic 
limit U(q) = provides the known moments of the harmonic-oscillator ground state. The 
C n = 2~ n [2] are then fixed, and we have 

r a, n _ (n-a)\a\ f U^A {2a ~ n)/i 

^o,o - 2n((n _ a )/ 2 )!(a/2)! V ma; 2 J [ ] 

for even a and n. 

As recalled here, the derivation of ( 120]) with its precise coefficient requires additional 
assumptions about the initial values of the moments, or about the kind of states whose 
evolution is considered. General effective equations are not unique owing to the depen- 
dence on classes of states described by them. The usual low-energy effective action (JXJ) 
is parameter-free only because it refers to a specific regime of states near the interacting 
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vacuum, as indicated by the qualifier "low-energy." The underlying conditions are tanta- 
mount to requiring the moments to agree with those of the harmonic-oscillator ground state 
when U(q) = 0, and indeed canonical effective equations with this choice are equivalent to 
equations of motion that follow from the low-energy effective action [2]. The solutions for 
moments then amount to expanding around the adiabatic vacuum state of the anharmonic 
system. Since our results build on (120]) . we will be dealing with the same states, but to 
higher orders in the semiclassical and adiabatic expansions. 



3.1.1 Solutions at zeroth and first adiabatic order 

Using zeroth-order solutions in Gq'q = {Gq'™, Hq}, the equation of motion at first order in 
A and zeroth order in y/h, we have, for odd a or n, 

G- a,n r> s~ia—l,n , l \ /~ia+l,n U \Q)^ s-ya—l,n /r>i \ 

0,0 — — —auj& 1 ' + [n — ajuG^ ' G 01 ' (21) 



But this equation is identical to ( 1181) . so we have the same solution, namely zero, for odd 
n. This pattern continues to all orders in the adiabatic approximation: 

= for odd n. (22) 

For even n, however, solutions change with progressing adiabatic order. We can still 
use (121]) for a odd and n even, describing solutions Gq'i with even a and n. Again, the 
equation is identical to (TT8|) . solved by ( TT9]) with new coefficients C' n instead of C n . If 
we match with the harmonic-oscillator ground state, we must require its values for the 
moments G%' q + G^'J to the present order, so that C n + C' n = 2~ n . Only this combination 
of C n and C' n appears in equations of motion to first adiabatic order, and therefore it is 
not necessary (nor possible) to determine both coefficients independently. Since we already 
fixed C n , we keep this value as well as Gq'q as in (120]) . This choice implies C' n = and 
therefore 

Gqi = for even a and n . (23) 

For moments of odd a and even n, finally, the solution of f[2~T]) (with odd a inserted) 
is different. In this case, the time derivative of ( 1201) becomes the left-hand side of the 
first-order adiabatic equation of motion: 

in - a)\a\ 2a -n U"'(q)q ( U"(q) \ ( 2a ~ n " 4 )/ 4 



2 n {{n-a)/2)\{a/2)\ 4 mu 2 \ mu 2 

— qojG 1 ' + (n — o)wG i ' G x ' (24) 



moo 

We can solve this immediately for a = n, and then substitute the result in to solve for the 
case a = n — 2, and so on. The general solution is 

U"'(q)q ( U"(q)\ ( - 2a ~ n ~^/ 4 
Gq'" = C ain — I 1 H — ) for odd a and even n (25) 



moj 3 \ mcu 2 J 
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where the C a , n are dimensionless prefactors. In particular, while the zeroth adiabatic order 
does not allow quantum correlations, they may appear at first adiabatic order, for instance 
by G 1 ' 2 ^ 0. 

In the harmonic limit, (125j) vanishes identically, and therefore the C a>n are not restricted 
by the requirement of perturbing around the harmonic ground state. Instead, the C a ^ n 
are fully determined by the adiabatic equations. For a = n in fl24|) . we have C n _i in = 
_2-("+ 2 ) n !/(n/2)!. Plugging ((25]) into fl2|) we find 

n — a (n — a)\{a — 1)! 

Ca ^ n = — Ca+1 > n ~ 2^{{n-a)/2)\{a/2)\ [2a ~ n) ' (26) 

a recurrence relation solved by the general expression 

(n — a)\(a — i)\ . 

Ca - 1 ' n = ~2«+2((n-a)/2)!(a/2)! (2a_nJ (27) 



(n-a-2)/2 

-n-2 



E 

b=0 



nn — (a + 2c) 
a + 2c 

c=0 



(n-a^-l)! 
((n-<)/2)!«/2)! 



(2a' b - n) 



for even a, where a' b = a + 2(6 + 1). 

To summarize, at first adiabatic order only moments with odd a and even n change, 
depending on the time derivative of q in ( 125]) . 

3.1.2 Second Adiabatic Order 

At second order in A and zeroth order in h, the equation of motion is 

GZ = -auj (l + G a ^ n + (n - a)uG a ^ n . (28) 

For odd n the solution is zero by (122|) . and even a and n in the equation of motion leads 
to Gq'2 of the form (TL9I) for odd a and even n, again with new coefficients C n . For these 
moments to vanish in the harmonic limit, we have Gq^ = for odd a. 

For odd a and even n, the left-hand side is given by the time derivative of equation 
. Substituting in this result and rearranging slightly, we have 



r °< 2 " n -a V mw 2 I °< 2 



+ 7 T i (t^(g)ff + EH?)? 2 } (i + ^l 



(rt — ajma; 4 \ \ mar / 



(2o-n-10)/4 



,2a -n- 6 /„ E/"( )\ 
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The form of this equation suggests an ansatz 



-\{U"\ q )q + U""{q)e) (l + ^f) 
1 ma; 4 \ mar J 



(2a-n-8)/4 N 



-Ba,n -r-^(U"\ q )qf 



^ , U"(q)\ 
4m 2 uj 6V ^ mco 2 J 

n/2\(n\- 1 U"(q)^ a/2 



(2a-n-12)/4 N 



a/2 J \a J \ mui 



,2 



Gjg (30) 



where A ajn and 5 a n are dimensionless coefficients determined by a and n, and the final 
term is motivated by the solution at zeroth order in A, generated by the first term in ( l2T?j) . 
Using this ansatz in ( |29|) . we find recursion relations for A an and B an : 

A a +l,n — : 1 A.a-\,n (31) 

n — a n — a 
C a , ra (2a - n - 6) o 
n — a n — a 

for odd a. Consistency of fl3"U|) requires that A n = B 0n = 0, which gives A,n = ^zy and 
B 2t n = — ^rf Ci,n- The recursion relations, with these initial values, give 

(a-3)/2 / 6 \ 

^ = ^+ E II T 9 J T^T^m ( 33 ) 

n_a to V=o ~ ( a - 2c )/ «- (a -2(6 + 1)) 

for odd a > 3, where the C a ^ n are given in ff26|) . The expression for the S a;n is similar, the 
only change being that C a , n is replaced by (2a — n — 6)C ajn in each term. 

To fully determine the moments, which are thus far given in terms of Gq' r 2 \ we need a 
condition from third order. Since we are still at 0(h°), the third-order equation of motion 
has the same form, so the condition is given by [I] 

even a v ' ' v ' 

for even n. From fj30|) we see that this is a complicated differential equation for Gq' 2 . Given 
(1301) . this condition suggests an ansatz for Gq' 2 of the form 

gss = ^^on?)? + t/ /w (g)g 2 ) f i + ^y + ^-^(ch^) 2 f i + 

ma; 4 \ mw 1 / 4m z u b \ mcj z J 

(35) 

where r, s, A^, and B' n are some undetermined constants. (Given a differential equation 
( )34|) . we expect one free parameter, which would be multiplying all Gq 2 . However, since 
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we are also solving the recurrence relation (!30|) . whose first two terms are fixed, imposing 
the consistency condition (|34|) will not leave any free parameters.) 

Substituting ( 1301) into ( )34l) with (j35l) . differentiating, and moving things around a bit, 
we have that 



moo 



-T 



n/2 
a/2 



.4, 



d (£/ w g + U ""q 2 > 

dT 



Am 2 u 6 



a/2 



d(U"'q + U ""q^ 

dT 



X r+ 2 + r 



2a - n - 8 

4 ma; 2 

a\ fT'g 



(U"'q + f/"V) 



2/ muj z 



(U"'q + U""q 2 )X r+ %- 1 



n/2 
a/2 



5, 



2U'"q (U"'q + U""q 2 ) + 2a ~ n ~ 12 U "'\ {U »>qf X ^ 



, (n/2 



n \a/2 



4 

2U"'q (U"'q + U""q 2 ) X s+ % + ( s + ^)^ {U'"qf X 
v ' V 2/ mar 



mar 



(36) 



must vanish, where X = 1 + U"/mu 2 . By inspection, we see that for r = — (n + 8)/4 and 
s = — (n + 12) /4, terms involving the same expressions with q and its derivatives also have 
the same power of X. This leaves only the numerical coefficients to be fixed. Only the 
first terms in the first two lines are proportional to X^ n ~ & ^ 4 d(U"'q + U""q 2 )/dt; generically, 
these terms must add to zero separately, which allows us to solve for 



A' 



V ( n/2 )A 
V ( n / 2 \ 2 ( n \~ 1 

Z-~ia \a/2J \aJ 



(37) 



Similarly, only the last terms in the last two lines are proportional to (U"'q) 3 X^ n ~ 16 ^ 4: , the 
remaining terms being proportional to U"'q (U"'q + U""q 2 ) X^ n ~ 12 ^ 4: , so we can solve for 
B' in similar fashion: 



BL 



E a C; 2 2 Rn(2a-™-12) 



(3* 



To confirm that these expressions are indeed valid, we need to check that the remaining 
terms in ( l36l) vanish. Substituting in our expressions for A' n and B' n , factoring out common 

quantities, and making use of the fact that Y2 a (a/2) 1 (^ a — n ) = 0, we find that 



remamm 



ining terms oc ^ (a/2) ^° ~ ^ ^ Aa ' n + Ba ' n ^ ' 



(39) 



We have checked that this expression vanishes for n = 2,4, 6, confirming the solution 



G 



O.n 
0.2 



-^- 4 (U"'(q)q + U""(q)e 



U"(q) 



-(n+8)/4 



B'n 

4m 2 w 6 



(U"'(qW 



U"(q) 
muj 2 



raw 

-(ra+12)/4 



(40) 
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at least to these orders, where A' n and B' n are given above. This expression also reduces to 
the solution for Gq 2 given in [I], with the correct coefficients A' 2 = 1/16 and B' 2 = —5/16, 
which can be checked using the earlier expressions for A a:Tl , B a>n , and G Qjn . The solution 
for Go'2 does not modify the ground-state condition for C n in ( 1T91) because it automatically 
vanishes in the harmonic limit U(q) = 0. Instead, the coefficients A' n and B' n are completely 
fixed without a choice of state. 

3.2 Adiabatic Approximation at O(Vh) in the Moments 

Starting with the first order in \/h, we need to consider different adiabatic orders in sepa- 
ration. 



3.2.1 Zeroth Adiabatic Order 

At leading order in the adiabatic approximation, {Gg'g , Hq} = 0. At O(Vh), that is, e = 1, 
this gives (from Eq. ffT4"j) ) 



n a-\,n . / \ n a+l,n U"(q)a na -l,n . U"'(q)a ^-,0,2^(0-1^-1 

—auLr 10 + [n - a)ULr 10 l, 10 + — rr^^o.o^o.o 

mu 2{muj) 6 ' z 



U"'(q)a ( „ a _i, n+ i (a - l)(a- 2) a _ 3iW _ 3 \ 

^0,0 To ^0,0 ] • l 41 J 



2(mu) 3 /2 V ' 12 

The solutions at O(h ) are given by ( 1201) for even a and n, and Gq'q = for odd a or n. 
We can use this result in ()4ip to obtain solutions for the G^'q. For even a or n in (|4ip . all 
the Go,o terms vanish, and the equation is identical to f fTBj) . giving the same solution for 
the moments involved: 

Gl$ = for odd a (42) 
(n-a)\a\ ( _ U" {q)\ {2a - n),A 



Gi'n = C'— . . ... — — - 1 H — for even a and n (43) 

1,0 n {{n-a)/2)\(a/2)\ \ mu 2 ) v ; 

As before at first adiabatic order, we implement the ground-state condition by requiring 
Gqo+V%G"'q to agree with the known harmonic values when U(q) = 0; thus, Cn + y/hC" = 
2~ n (G" is not dimensionless) . Again keeping Gq'q unchanged compared to fl20|) . we have 
G" = and therefore GI'q = for even a and n. 

For odd a and n in (j4"Tl) . the Go,o terms do not vanish. Substituting fl20|) into fj4Tl) and 
simplifying the resulting expression, we have 



n / \^ya+l,n ( -. . U (o) \ ^o— l,n 

= (n - « G li0 - a 1 + r 

£Hg)a 4a-3n-l /a\ / n-a + l \ / t/"( g ) \ (2a ~ n " 3)/4 
+ m 3/2 w 5/2 12?r 1 v27 V 2 J \ mu; 2 J ' 1 1 
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For the case a = n, this gives 

C n-i, n _ U"\ q )a n-l n fl\ f U»{q)\W 

We can plug this solution into the a = n — 2 equation to solve for G™Q 3 ' n , and so on. In 
general, 

Qa-i,n _ n — a ( U"(q) \ ^, a +i,n 
'° a \ moo 2 J '° 

U"'(q) 4a-3n-l ,a\ ( n-a + l \ ( U"(q) \ ( 2 — 7 V 4 

From ( l4"5j) . we see that the two terms have the same power in 1 + U"(q)/moo 2 , which 
allows us to write the solution as 

U"'(a) ( u"( Q \\(^-n-5)/4: 

Gt'o = Da,n 1 + fOT eVen a aIld ° dd n ' ( 4? ) 

m A i l bj b i l \ moo 2 J 

where 

f-l) 6 r (z) f ( 1\ 

= k(T=li(("- 1 ' 6! ^ + ("- 86 - 1 ' r ( 6+ 2) 

- - 8(6 - c - 1) - l)r ^6 - c - (-6) 

if n > 5 and b > 2, and 

12vr(n-"2) 



c+1 



(4? 



£>a,n - <| JSnJt T j|j f if n > 3^ 6 = 1 ^ 

is a dimensionless prefactor that depends on a and n. In the above expression, b = 
(n — a — l)/2 and (x) n = x(x + + n — 1) is the Pochhammer symbol. Comparing to 

(!20|) . we see that for odd n and even a, 

G 5?«( 1 + ^) '^T 1 - (so) 

The additional dimensionful factor of U'" /m 3 ^ 2 oo 5 ^ 3 in (j4"T|) provides the correct dimension 

of h" 1 ' 2 . 
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3.2.2 First Adiabatic Order 



At first order in the adiabatic approximation and at 0{y/%), the equation of motion is 

Aa,n / \ < ^o+l,n /, . U"(q)\ a _ ln U"'(q)a f „Q2 n a-\,n-\ 

G i,o = [n- apG^ -aw 1 + G lt + — G 00 G 01 

_,_x-y0,2x-ya-l,n-l _ ^a-l.n+l , ( Q ~ 1)( Q ~ ^) n a-Z,n-1, 
+ ( - r 0,l ( - r 0,0 ^0,1 i ^ o. 1 

If a is odd, we see from ( H2l) that the left-hand side is zero. If n is even, we see from ( [22]) 
that the last four terms on the right-hand side vanish. Let us consider the simplest case 
first: odd a and even n. In this case the equation is identical to f|T8|) . giving the same 
solution f|T9|) for the moments involved, with a new (and dimensionfull) C n . Recalling that 
the 0{h 1/2 \°) and 0{h°\) solutions (g3} and (p} were also the same, we can write 

s-ia,n r^a^ri s-ia,n s-ia,n 
b l,l K ^0,1 ^ ^1,0 ^ °0,0 

Now requiring Gqq +Gq'" + v^(G°'q to agree with the known ground-state moments 

in the harmonic limit and keeping Gq'q as in fl20l) . we have 

Gl'i = Gq'i = GI'q = for even a and n . (52) 

In these equalities we have made use of another relation. It turns out that the simpli- 
fication of the equation of motion observed here is by itself not sufficient to guarantee the 
same solution, because the equation does not fully determine all the moments. As in the 
zeroth-order approximation in both y/h and A, a constraint ( ]34j) is needed from the next 
adiabatic order in order to determine G 0,n [I] . (More specifically, the constraint shows that 
C n does not depend on q.) Here, the second-order adiabatic equation of motion at 0(\/h) 
is 

W,i -{n- a)ujG 12 -awli + ——^ j G 12 + 2 ^^ 3/2 ^ ,o G o,2 

, ^,0,2^0-1,71-1 ^-yO^^-ya-l.n-l _ ^-ya-l,n+l . ( a ~ 1)( Q ~ 2) ^-3^-3 
+ Lt 01 Lt 01 -+- <-t 2 Lt Lt 2 -+- ^0,2 

But n is still even, so the last five terms vanish, again due to ( I22p . The right-hand side 
once again is the same for all cases considered for now, so the same condition on the left- 
hand side follows as in [I], and (152]) is indeed correct. We note, however, that not all these 
equalities between the moments will be valid at second order in A because G^i, appearing 
at the left-hand side of the equation of motion (I5U|) . will no longer be zero for odd a and 
even n. 

Now let us consider the case where a and n in ( 151]) are both even. The extra terms 
on the right-hand side vanish, but the left-hand side is given by the time derivative of 
G"'q . The resulting equation, which describes moments with odd a and even n, is identical 
to ( 124]) . the corresponding equation at 0(h°) — except that the coefficient 2~ n is the C n 
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belonging to G^'q, with G\'q = for even a and n. In this case the equation of motion 
fully determines all the moments; no constraint from the next order references to harmonic 
states are needed. Consequently, the solution is 

Gi'i = for odd a and even n. (54) 

For odd n, the solutions are not the same as the 0(h°) solutions ([22]), which are zero. 
In the case where a and n are both odd in ( 151]) . the extra terms remain, but the left-hand 
side is zero. If we compare ([5"T]) to (14 ip and recall ( ]52j) . we see that the additional terms 
are identical, except that the G°' 2 G a ~ l,n ~ 1 term in ( I4T1) effectively occurs twice in ( IBTT) . 
This has the effect of replacing the factor of 4a — 3n — 1 in ( 14*41 with 4a — 3n + 2, so the 
solution is the same up to a change in the prefactor in ( 1471) . 

The final case, where a is even and n is odd in ([ST]) is more difficult. The left-hand 
side is nonzero and given by the time derivative of ( I4"T|) . and the extra terms on the right- 
hand side involve the expression given in ([25]) . Here we will not attempt to find a general 
solution. 



3.3 Adiabatic Approximation for n = 2-moments at 0(h) 

At second order in \/h, the equations again get more complicated. We will restrict ourselves 
here to deriving a relation for solutions of second-order moments (n = 2) as needed for the 
leading correction in the equation of motion (jTTj) . Our considerations in this section only 
illustrate the procedure but do not provide complete solutions. 

The n = 2-moments, which are the same at 0(h°) and 0(a 1/2 ), are different at 0{h). 
The 0(h) equation of motion for the moments is given by [2] 



r 2 = — awG 2 + [n — cl)ojG 2 G 2 



7TILU 



U" (q)a i n ^i n a-l,n-\ , n 0,2 n a-l,n-l\ 



2{muf/ 2 



U"'{q)a ( G a-i,n+i _ (a - l)(a - 2) ga „ 3iW _ 3 



2(mu) 3 / 2 V 12 

U""{q)a „Q3„ a _ hn _ x U""(q)a ( 1+2 (a-l)(a-2)^. 



+ 3!(^F Go G ° " 6(^)* V 4 G ° J (55) 

At zeroth order in the adiabatic approximation, the left-hand side is zero. For a = 0, we 
find that G20 = 0; as &t 0{h 1 ^ 2 ). In particular, G^q = 0. The a = 2-equation gives no new 
information, but confirms that G\'q = 0, as we found in Sec. 13.2.11 The a = 1-equation is 

n _,, r 2,2 U"'(Q) ^0,3 ^0,4 ,, fi x 

° " UG ^° ~ u [ 1 + ^JT ) G 2,o - 2(mw) 3/2 G i,o - "^p G o,o • (56) 
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The coefficients G^'q and G^g are gi ven by equations (|47|) and (1201) . respectively, so we can 



express G 2 . i n terms of G^'l'- 



U"(q)\ ^0,2 {U'"{q)f / £/"(?) V 2 , f/""(g) A £/"(<?) 



^"V 1 m^'J^o V 24m 3^5 ^ mc ^J + 8rr ^( 1 ^) • (57) 
In order to find G^'q, we would have to go to first order in A. 



4 Equation of motion for the oscillator up to O 
and the fourth adiabatic order 

It is evident from (I4UI) compared to the previous equations that higher time derivatives 
appear at higher orders of the adiabatic expansion. This pattern continues at higher 
adiabatic orders, even going beyond the second derivative order of the classical equations. 
Genuine higher- derivative equations are then obtained when solutions for moments are 
inserted in ( IT2l) . Not surprisingly, it becomes more and more complicated to find explicit 
solutions to higher orders, valid for generic a and n. Still, individual moments, such as G ' 2 
as needed for the first corrections in ( 1TTT) . can be computed more easily because specific 
numbers take the place of coefficients such as A' n and B' n subject to complicated recurrence 
relations. 

4.1 Higher-derivative equation of motion 

From the preceding section it is clear how such calculations are organized, and it suffices 
here to quote the results needed for the leading corrections in the equations of motion of 
the anharmonic oscillator, given as before by 

q = m~ 1 p 

00 1 / 1 \ n/2 

p = -mu 2 q-U\q)-Y ~, \— T U (n+1 \q)G°' n . (58) 

^— ' n! V mcun ) 

n=2 x ' 

Taking the time derivative of the q equation, we may write, correct up to 0(h 3/2 ) 111 
quantum corrections, 

q=-UJ 2 q-U\q)/ m --^ 



2m? 



ijj 



showing which moments and orders we need. (Moments of orders higher than n = 2 would 
be required at the next order, 0(h 2 ).) Here we have already used the fact that Gq^ = 
for any value of i, according to (|2"2~|) . 

In order to evaluate (I58"|) completely, we need to compute G ' 2 to orders 0(h°) and 
0(h 1 ^ 2 ) in the semiclassical expansion. The previous section contains results up to the 
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second adiabatic order, at which the right-hand side of (|59|) would be of second order in 
time derivatives. In particular, we have 

from (1201) . while 

G?J = G?g = GS3 = 0. (61) 
At second adiabatic order, we have 



,0,2 = U"'(q)q + U""(q)q 2 L + U"(q) \ ^ _ h{U"\q)qf / | U"{q)\ ~ 7/2 
16mu; 4 V mu 2 J QAm 2 u e \ mcu 2 



while 

G° { 2 = . (63) 

Higher adiabatic orders cannot be obtained from the previous general formulas for 
G a ' n , but the relevant contributions to the specific moment G ' 2 can be computed with 
the previous methods. Up to fourth adiabatic order, we have Gq 3 = (which holds for 
all moments of even a and n) by using ( 1281) at higher adiabatic order and observing that 
Go'2 = for odd a and even n. Moreover, G^ 2 = and G a {\ = vanish. The final moment 

2 

we need, G ' 4 is more difficult to derive, and we just sketch the procedure we followed. The 
fourth-order equations can be manipulated to show a relation of the form 

Gli = XG° > 2 + ® (64) 

where X = 1 + U"(q)/mu 2 and 



+ 

+X -9/2 



( 15 



{ 64^7 J [U "' m t U '" (q) * + U '"" m * + 3tr 'M 



_ xll/2 / 315 \ 
^SUm^u 1 ! J 



[U'"(q)q] 4 (65) 



Then using the consistency equation (13~4|) from [2], we get 



2XG° 'l + ( ^f] G° 'l + ±0 = . (66) 
18 



Choosing G° ' 2 4 = AX~ 7/2 + BX~ 9 / 2 + CX" 11 / 2 + DX^ 2 as an ansatz, we solve for these 
functions A, B, C and D. However, we have five equations (corresponding to the five 
different powers of X in ( l66il ) with four unknowns. This generates another consistency 
equation which turns out to be satisfied by the coefficients in 



G o i2 U'"(q)"q + AU""(q)qq + 3U""(q)q 2 + W""(q)q 2 q + U"""(q)q 4 / U"(q) ^ 

6Amu 6 \ muj 2 

2l{U""{q)q 2 + U"'{q)q) 2 
256m 2 w 8 

7(U"'(q)q)(U"'(q)q + 3U""qq + U""'(q)q 3 ) \ f U"(q) X ~ 9/2 
64m 2 uj 8 J \ mco 2 

231(U"'(q)q) 2 (U"'(q)q + U""{q)q 2 ) / £/"(g) N _11/2 



512m 3 o; 10 \ mw 2 

1155([/ w (g)g) 4 A , U"(q) \ " 13/2 ffi7 , 
4096m 4 w 12 V m^ 2 / ' 

We may now rewrite the equation of motion ( 158|) as: 

g = —cu 2 q — U'(q)/m (68) 
£/"'(?) [/(5,g) + /i(?,g)?+ / 2 (g)<f + /s(?,?)9 + /4(g)'g"] 



1m 2 uj 



where 

/(g,g) 



t^M \ ~ 1/2 + ^ /w (g)g 2 A + u%qY\ ~ 5/2 _ 5(fA'( g ))V A + tr(g) 

2 \ mw 2 / 16mcj 4 \ mu 2 J 64m 2 w 6 \ mu 2 

U"""(q)q 4 f i + U"(q) \~ 7/2 + 21(U"" \q)) 2 q 4 / , £A(g)y 9/2 



-7/2 



64mw 6 \ mw 2 / 256m 2 oj 8 \ mw 2 

7U""'(q)U"'(q)q 4 + tA(g) \~ 9/2 _ 231fA"(g)(tA'(g)) 2 g 4 A , A"(g)^" 11/2 



64m 2 oj 8 \ mu 2 / 512m 3 co 10 \ mcu 2 

1155(U"'(q)) 4 q 4 A t/"(g) \ " 13/2 



4096m 4 c- 12 * — 2 



(69) 



^ C/ W (g) A , t/"(g)V 5/2 3U"">(q)q 2 A , C/"(g)^ 



16mo; 4 \ mw 2 / 32mu; 6 \ my 2 

1 63^ w (g)[/ w (g)g 2 A + ^(g) A ~ 9/2 _ 23i(^( g )) 3 g 2 / + u"( q ) \ ^ n/ ( 2 7Q) 

128m 2 u; 8 \ mu 2 ) 512m 3 co 10 \ mw 2 / 
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/a(g) 



3U""(q) 
64mw 6 

U""{q)q 
16mu 6 

U"'(q) 
64mu 6 



1 + 



U"{q) 
mu) 2 

U"(q) 
mui 2 

U"(q) 
mu 2 



-7/2 



-7/2 



2l{U"'{q)) 2 
256m 2 ui 8 

W'{q)?q 

64m 2 w 8 



U"(q) 
mu> 2 

U"(q) 
mu 2 



-9/2 



-9/2 



-7/2 



(71) 
(72) 
(73) 



Once these coefficients are inserted in (I58I) we have the equation of motion for q correct up 
to the fourth adiabatic order, for quantum corrections up to h 3/2 . From flSZD, it is now clear 
that higher time derivatives result, up to fourth order with the present approximation. The 
equations for moments shown in this section demonstrate that it is the adiabatic order, 
rather than the semiclassical expansion or the order of moments, that determines the order 
of derivatives. Although back-reaction of moments on expectation values is responsible for 
higher time derivatives, the new degrees of freedom that higher-derivative equations would 
imply if used at face value, are not identical to the moments as true quantum degrees of 
freedom. 



4.2 Uncertainty Relation and Zero-Point Energy 

Having derived solutions for second-order moments, we need to make sure that they obey 
the uncertainty relation so that they can correspond to a state. This requires us to compute 
not only the moment G ' 2 , as used in corrected equations of motion, but also G 1,2 and 
G 2,2 . Since expressions with these moments get more lengthy, we restrict the orders to 
(0, 0) + (0, 1) + (1, 0) + (1, 1) + (0, 2). (Especially G 1 ' 2 is more difficult to obtain to higher 
orders.) We then have 



G o,2 = 1 L , U"(q)Y 1/2 , U"'(q)q + U""(q)q 2 L | U"{q)^ 
2 \ mu 2 J 16mo; 4 \ mu 2 

5{u"'{q)q) 2 A , u"( q )Y 7/2 



64m 2 uj 6 \ mu 2 
G 2,2 = l( 1+ U"(q)\ 1/2 U"'{q)q + U""{q)q 2 (, , U"(q)y 3/2 



(74) 



2 y mu 2 J 16mu; 4 V mu 2 

7(U"'(q)q) 2 / , ^7"(g)^- 5/2 



+ Ha 2 6 1 + (75) 
U"'{q)q(. U"(q)Y 3/2 



Qmuj 6 \ mou z J 

To check the uncertainty relation G°' 2 G 2 ' 2 — (G 1 ' 2 ) 2 > 1/4, it is useful to write G ' 2 = 
l x -i/2 + Y w hile G 2 ' 2 = (iX"V2 - y + 1 S0 ! ^" 7/2 ) X, where X = 1 + ^ and 

y = lT{q)q + U""{q)q 2 , U"(q)y 5/2 5(U"'(q)q) 2 , U"(q)y 7/2 



lQmu 4 \ mu 2 I 64m 2 cu 6 \ mu 2 
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With these definitions, the required expression is of the form 



0,2.^2,2 



G^ 2 G 



(G 



1,2\2 



X 



ix^-Y^ + ax-^ + Yy 



32m 2 w 6 



64m 2 w 6 



1 /4 - XY 2 + 



32m 2 u 6 



X^ 2 Y. 



(77) 



In the limit U(q) — > 0, we have X = 1 and Y — 0, yielding a value of 1/4 for the 
above expression and exact saturation of the uncertainty relation, in accordance with our 
assumption of the Gaussian ground state in the harmonic limit. With anharmonicity, it 
is not obvious to see the sign of G°' 2 G 2,2 — (G 1 ' 2 ) 2 — 1/4. It is, however, clear that we 
need Y > 0. The adiabaticity condition for the anharmonic ground state is therefore not 
guaranteed to be valid automatically, but with our equations the uncertainty relation can 
easily be monitored when numerical solutions are analyzed. 

Finally, having determined the moments G 0,2 as well as G 2,2 to some orders, we can 
compute anharmonic corrections to the zero-point energy Z = ^hoj(G 0,2 + G 2 ' 2 ) =: ^huZ'. 
With (I74|) and (I75p we have, now valid up to the order (1,3), 



Z 1 



G 



0.2 



G 



2,2 



U"(q) 



-1/2 



U"'{q)q + U""{q)q 2 



5(U"'(q)qf 
G4m 2 u & 



16mo; 4 
U"(q)^- 7/2 



in 



1/2 



U"'(q)q + U""(q)q 2 



64m 2 w 6 



U"(q) 



16mcu 4 

-5/2 



U"{q) 



muj 1 



U"{q) 



-X^ 2 (l + X) + Y(l -X)+ (^ W (g)g) 2 x -5/2 
2 v ; v ; 32m 2 w 6 



-5/2 



-3/2 



(78) 



5 Conclusions 

We have shown explicitly how higher time derivatives arise in effective equations of canon- 
ical quantum systems. Physically, quantum back-reaction by moments of a state on ex- 
pectation values implies correction terms, which in an adiabatic approximation (combined 
with a semiclassical one) can be solved for in terms of higher time derivatives. The quan- 
tum degrees of freedom are then removed from direct view, analogous to integrating out 
variables in a path integral. Although the quantum equations are local, given by differ- 
ential equations of finite order, the effective system in which infinitely many degrees of 
freedom have been expressed by higher time derivatives becomes non-local. 
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5.1 General properties 

Our analysis has revealed several general properties of the expansions considered. 

First, some coefficients of moments vanish at all orders in one of the expansions, ir- 
respective of the harmonic state perturbed around. For instance, as already mentioned, 
Go'" = for odd n. This observation simplifies some computations of moments relevant 
for corrected equations of motion. The interesting moments that appear in quantum cor- 
rections, however, satisfy equations that become progressively more involved as the orders 
increase. The methods employed here can be extended to higher orders, but a general 
solution seems out of reach. 

Secondly, while coefficients of explicit solutions for moments are usually complicated, 
the relationship between the adiabatic and derivative order is evident from the general form 
of equations (|T6|) . Going to higher adiabatic orders requires the use of time derivatives of 
lower-order coefficients, starting with zeroth adiabatic order in which coefficients depend 
on q. In this way, higher and higher time derivatives enter solutions as seen explicitly in 
the preceding section. 

Thirdly, as in our explicit examples ffl9|) . (125]) and fl30|) . odd adiabatic orders mainly 
change moments G a,n with odd a (correlations), while even adiabatic orders lead to cor- 
rections in moments with even a. 

5.2 Higher time derivatives in quantum cosmology 

Our results apply to all quantum systems, but are especially relevant for quantum gravity 
and cosmology. In these settings, canonical techniques are often crucial or at least applied 
widely, and it has remained unclear if and how higher time derivatives should result. 
The present article shows this unambiguously and provides a systematic procedure for 
their computation. As an example for the importance of higher time derivatives, we may 
look at loop quantum cosmology [TUl E]. In this setting, one would generically expect 
effective equations with higher time derivatives, as a result of higher- curvature corrections. 
(We are dealing here with quantizations of finite-dimensional systems, and must therefore 
leave unaddressed the intriguing possibility that strong quantum corrections may lead to 
signature change from space-time to a 4-dimensional quantum version of Euclidean space 
[TTj [T2~| 113] . The question of how the "evolution" of states and quantum back-reaction 
can be formulated to compute quantum effects in timeless Euclidean space requires further 
detailed study. But quantum effective equations should still be non-local and require higher 
derivatives, including spatial ones.) 

However, higher derivative terms or the more general moment- dependent corrections 
as in (TIT]) have not yet been computed in all cases; adiabatic regimes may even be non- 
existing. (After all, non-perturbative quantum gravity may lack a ground state with slowly 
evolving moments.) Instead, for a first impression of implications of quantum effects one 
often ignored quantum back-reaction altogether, considering only quantum-geometry cor- 
rections in homogeneous models which are easier to implement by simple modifications 
of classical equations. A prominent one is the so-called holonomy modification, by which 
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the Hubble term "H 2 = (a/ a) 2 in the Friedmann equation is replaced by a periodic func- 
tion such as sin(£'H) 2 /£ 2 with a length parameter £ (which could be the Planck length); 
see e.g. [T4l IT5] . This modification is motivated by a property in the full theory of loop 
quantum gravity [Pol ITT] , according to which only holonomies but not the gravitational 
connection can be represented as operators. A heuristic interpretation often encountered, 
states that higher-order terms in an expansion of sm(£'H) 2 /£ 2 = H 2 (l — \£ 2 1-L 2 + 0((fH) 4 )) 
are related to higher-curvature terms. However, this interpretation overlooks the fact that 
higher-curvature terms also provide higher time derivatives, which are not included in most 
studies (and have not yet been computed). But if the higher-derivative part of curvature 
corrections is ignored, one cannot consider isolated higher powers of spatial curvature com- 
ponents as a reliable expansion. Generically, there is no reason to assume that a term of 
"H 2 is more important than, say, both of which contribute to the space-time curvature 
scalar. Expansions become inconsistent when only one type of terms is kept and, perhaps 
even more damningly, general covariance is put in jeopardy. 

With holonomy corrections, the whole series expansion of sm 2 {£%) / £ 2 by £H is used, but 
not a single higher time derivative. Such an approximation cannot be consistent unless one 
can show that there are no higher time derivatives whatsoever. There is indeed a harmonic 
system in loop quantum cosmology free of quantum back-reaction [6], given by a spatially 
flat isotropic model with a free, massless scalar. In this model, holonomy corrections 
correctly describe quantum evolution. If one departs from this model just slightly, when 
matter remains kinetic dominated and anisotropics and inhomogeneity are small, quantum 
back-reaction does arise [T8l IT9] but may be assumed weak; holonomy corrections can still 
be reliable. In all other cases, however, correct physical conclusions can be drawn only 
when all quantum corrections have been computed. For instance, when there is a phase of 
slow-roll inflation, the potential dominates and one must expect effective equations based 
on the assumption of kinetic domination to break down. The presence of higher-derivative 
terms means that holonomy corrections on their own are not reliable, but it also suggests 
interesting relationships with early-universe models based on higher-curvature or non-local 
derivative terms, especially regarding the singularity issue [201 12H 122] . There is still much 
work to be done to apply the complete effective methods to quantum cosmology. At the 
very least, the present article serves to clarify the role of higher-derivative corrections. 

5.3 The Role of Adiabaticity 

The use of the adiabatic approximation as the key ingredient to arrive at higher time 
derivatives implies that not all quantum regimes may be amenable to a higher-derivative 
description. The validity of the adiabatic approximation is an assumption, which can be 
tested self-consistently but need not always be valid. We had to use equations of different 
orders to determine all coefficients, and depending on the given quantum dynamics, not 
all these equations may be mutually consistent, for instance if anharmonic constructions 
are attempted for fully squeezed, correlated harmonic coherent states, contradicting the 
condition Gq'q = derived here for odd a and even n. Moreover, the final solutions 
may be mathematically consistent, but could violate the uncertainty relation. If they do, 
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there would be no state corresponding to the moment solutions. (See [23] for examples 
in quantum cosmology where these problems occur.) Solutions for moments describe a 
dynamical state, and not all states may allow moments to evolve adiabatically. Here, we 
have expanded around the ground state, which one can reasonably expect to evolve slowly. 
In other, more excited states, the condition may not be met. 

If the adiabatic approximation does not apply, there are still effective equations, and 
they can be expanded in a semiclassical approximation. However, it is then no longer 
possible to solve for the moments in terms of expectation values, and no higher-derivative 
effective equations exist. One would rather work with a higher- dimensional effective sys- 
tem, a dynamical system in which the expectation values together with all moments rele- 
vant to a given order in h are kept as in (II ip and ( fl2|) . The non-locality of the quantum 
system is then realized by the infinite number of moments if all orders are included. To any 
finite order in the moments, such a system can still be solved approximately, most often 
by numerical means, and evaluated for physical information; see e.g. [2U E] for examples 
including a rather large number of moments. Effective equations in terms of moments, 
explicitly exhibiting the true quantum degrees of freedom, are therefore more general than 
higher-derivative effective equations, but higher-derivative equations, if they exist, can 
show some features more directly. 
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